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We test emission models of circum- nuclear dust torii around quasars, at low and 
high redshifts, by using a large collection of photometric data for an unbiased sample 
of 120 optically-selected objects with millimetric and sub-millimetric fluxes, including 
new unpublished data. Under the assumption that the dust is heated by a point-like 
source with a power-law primary spectrum, as defined by the observed optical-UV 
continuum, we infer the basic model parameters, such as dust masses, temperature 
distributions and torus sizes, by numerically solving the radiative transfer equation in 
the dust distribution. In addition to the substantiated statistics, an essential improve- 
ment over previous analyses comes from the use of optical-UV data to constrain the 
primary illuminating continuum, which is needed to estimate dust temperatures and 
sizes. Dependences of the best-fit parameters on luminosity and redshift are studied 
and the contribution of dust in the host galaxy to the observed fluxes is briefly men- 
tioned. This analysis constrains the properties of the enriched interstellar medium in 
the galaxies hosting the quasars. The dust abundance does not display appreciable 
trends as a function of redshift, from z ~ 1 to almost 5, and shows that dust and 
metals are at least as, and often more, abundant at these early epochs than they are 
in local galactic counterparts. This evolutionary pattern is remarkably at variance 
with respect to what is expected for disk galaxies, like the Milky Way, slowly building 
metals during the whole Hubble time. It rather points in favour of a much more active 
phase of star-formation at early epochs, probably related to the formation of massive 
spheroidal galaxies. 

Key words: AGN - active galaxies: ISM, photometry - ISM: dust, extinction, con- 
tinuum - infrared: galaxies 



1 INTRODUCTION 

The presence of dust in high-redshift objects raised consid- 
erable interest in the recent years for mainly two reasons. 
First, there is a need to understand its effects on the visi- 
bility of the distant universe and possibly to correct for the 
biases this implies. Dust reddening of the light from distant 
galaxies biases all optical estimators of the star-formation 
rate and modifies the spectral energy distributions (SED) by 
mimicking an old stellar population. It has been suggested 
that the early bright star-formation phase predicted by some 
models for spheroidal galaxies could have been extinguished 
by du st produced by the fir st generations of mass ive stars 
(e.g. ( [Franceschini et al 199^ ), ( [Pe Zotti et al. 1996| )). Prop- 
erties of the dust distribution have then to be evaluated to 
correct for these various effects. 

In a more positive vein, observations of dust emission 
from high-redshift objects provide a vehicle of relevant in- 



formation on processes of stellar activity and metal yield in 
primeval objects, to be framed onto evolutionary models of 
galaxies. 

Indeed, not only evidence is accumulating that the most 
important episodes of star-formation in the local universe 
occurs in dusty environments (e.g. in the IRAS luminous 
galaxies), but also there are indications of an important 
role of dust in very distant objects. This is emphasized by 
studies of the hyperluminous IR galaxies at high-z (IRAS- 
F10214, IRAS15307, HRIO, Cloverleaf and objects detected 
in recent sub-mm surveys, Hughes et al. 1997; Hughes et 
al. 1998b), and by the fact that even optically ("drop-out") 
selected galaxies at z~2 to 4 (hence biased against dust ob- 
scuration) suffer dust extinction of 1-2 magnitudes at leas t 
i n the optical at such high redshifts (( Pettini et al. 1997 ), 
( Meurer 1997 )). Also, an integrated emission, in the form of 
an isotropic background, and corresponding to a very ener- 
getic high redshift phase, has been recently detected in the 
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sub-m m ((I'uget et al 199e|), (|Fixsen et al. 199^ ), ( |Schlege]| 
et al. I!p98[), ( [Hauser et al 1998| )). 

Ideal sites to look for a dust- and metal-enriched inter- 
stellar medium at high and very high redshifts are obviously 
provided by the environment of distant quasars, thanks to 
the enormous power available to illuminate and heat the 
surrounding medium. Indeed, the presence of an AGN could 
even mask the presence of a forming spheroid and compli- 
cate its investigation, but can also be exploited to identify 
target areas where to look for forming structures (galaxies, 
galaxy clusters and groups). 

Very high metal abundances are found in clouds pro- 
ducing the absorption line systems "associated" with QSOs, 
showing that the close environments of high-z QSOs were 
the site of quick and efficient star formation and metal pro- 
duction (Haehnelt & Rees 1993; Hamann & Ferland 1993; 
Franceschini and Gratton 1997). Quantitative analysis of 
emission and absorption lines would allow, in principle, to 
measure the metal abundances from the equivalent widths 
of atomic and ionic lines. However, these measurements are 
quite complex and uncertain, and particularly in the vicin- 
ity of the quasar the evaluation of the ionisation status is 
subject to systematic uncertainties. 

On the contrary, the physics ruling the generation of 
photons by thermal dust is much simpler and allows quite 
more robust and reliable evaluations of the dust mass. Of 
course, this implies an extensive sampling of that portion 
of the SED dominated by dust emission. Given the charac- 
teristic grain temperatures and the typical radiation fields 
around quasars, this is widely accepted to happen also in 
the FIR/mm domain of radio-quiet quasars (Sanders et al. 
1989; Danese et al. 1992; Hughe s et al., 1993, 1997; Rowan- 
Robinson 1995, (Lawrence 1997)). 



To this end, sensitive observations in the far-infrared to mil- 
limetric domain are essential. The improvement in sensi- 
tivity of modern bolometer detectors, combined with the 
favourable circumstance of the strong and negative K- 
correction implied by the steeply rising sub-mm spectra, 
have recently allowed to detect and characterize with sig- 
nificant detail the dust properties in distant objects. 

Many high-z sources, most of which associated with ac- 
tive nuclei, have been recently detected in the millimeter 
(both in the continuum and CO line emission; see Andreani 
et al. 1993; Chini and Kriigel 1994; Dunlop & Hughes 1994; 
Isaak et al. 1994; McMahon et al., 1994; Ivis on 1995; Ohta 
et al., 199 6; Omont et al., 1 996 a and 1996b; (Guilloteau et 
al 1997); ( [ivison et al 1998[ ); ( [Hughes et al 19971 )). Low-z 
AGN were observed by Chini et al. (1989a), Barvainis et al. 
(1992), Hughes et al. (1993, 1997) and the measured spectra 
nicely agree with those predicted by thermal emission from 
dust components at a variety of temperatures down to ~ 30 
K. 

In this paper we make use of a large sample of opti- 
cally selected quasars at low and high redshifts, with far-IR 
data by IRAS and all observed at mm or sub-mm wave- 
lengths. Since both detections and upper limits to the long- 
wavelength fiuxes are considered, the sample is to be con- 
sidered as an unbiased optically-selected one. 
We interpret the observed spectra in terms of a geometri- 
cal model in which the dust distribution is assumed to have 
simple axial symmetry, either confined to the quasar nucleus 
or covering a sizeable fraction of the hosting galaxy, and is 



illuminated by the central nuclear source. Obvious evidence 
in support of this is the dominance of the point-like nuclear 
source in the optical in both high- and low-redshift objects. 

Assuming that the radio-quiet quasars can be modeled 
as a homogeneous population and are representative of the 
AGN class as a whole, we build the entire infrared to m,m, 
spectrum via observations at fixed frequencies of targets at 
different redshifts. 

This analysis improves significantly on previous studies, 
as it exploits the whole spectral information from the optical 
to the mm,, the optical providing in particular an essential 
constraint on the intensity of the radiation field illuminating 
the dust. 

We adopt throughout the paper Ho = 50 Kra/s/Mpc 
and two values (go=0.5 and 0.1) for the cosmological decel- 
eration parameter. 



2 THE DATA 

The sample consists of 120 optically-selected quasars with 
sub-mm and m,m, observations and radio (either 6 or 20 cm) 
data. The basic selection criterion is the availability of long- 
wavelength (in addition to the IRAS) mm data and a quasar 
(rather than Seyfert galaxy) morphological classification. 

Tables 1 and 2 list the available flux data on the whole 
sample: optical magnitudes and near-IR fluxes (columns 
3-^8), IRAS detections or upper limits (columns 9-^12), sub- 
mm fluxes at 450, 850 and 1250 ^m (columns 13 and 14), 
radio fluxes (column 15) and absolute B magnitudes in col- 
umn 16. All fluxes are in units of mjy. Radio-loud objects 
are separately listed in table 2. 



2.1 Optical data 



Optica l magnitudes are taken f rom (gchneidcr et al 



(1991) ), ( [Schneider et al. (1992)| ), ([Storric-Lombardi ct al 
(1996)) and only for a few remaining objects from the 



( [Vcron fc Veron (1998) ) catalogue an d are corrected for in 
tergalactic absorption according to ( Madau (1995) ). The 
absolute B-magnitudes are taken from the same authors 
and for uniformity reported to go = 0.5. R-magnitudes 
are taken from the ESO/ST-ECF USNO A.l database 
(tittp://archive. eso.org/skycat/servers/usnoa). 



2.2 Near and Far-IR data 



NIR fluxes a re taken from ([Neugebauer et al. (1987) ), ( El 



etal. (1994)[ ), ( [Zitelh et al. (1993)[ ) and ( [Taylor et aL (1996)^ 
The far-IR data are taken from the IRAS PSC and FSC, 
and for most of the sources co-added fluxes were provided 
by IPAC based on the SCANPI (Scan Processing and In- 
tegration Tool) program. This procedure performs a one- 
dimensional coaddition of all the IRAS survey data on the 
source. The sensitivity is comparable to that achieved by the 
FSC (Faint Source Catalog) for point sources (see the IPAC 
manual for details). IRAS photometry for low-z an d some 



high-z objects hav e been taken from the literature ((Neuge- 



1992)) 



bauer et al. 1986), (Banders et al. 198!:), (Barvainis et a] 
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2.3 Millimetric data 

Mm observations are partly collected from the literature and 
partly obtained by the authors with the IRAM SOm at Pico 
Veleta. The low- redshift sample cons ist s mainly of PG radi o- 
quiet quasars (( |Chini et al 1989a|), (|Chini et al 1989b|)). 
High-z o n es from (( Chini and Kriigel 1994 ) and (|[saak et 
al. 1994|), ( pvison 1995p , ( [McMahon et al. 199^ , ( pmont et 
al. 1996b| )). 

New unpublished data (sources 0910+56, 0953+47, 
1158+46, 1247+34, 1548+46 and 2048+01) were taken by 
us with the MPIfR 7-channel bolometer (Kreysa et al., 1993) 
at the focus of the IRAM 30m antenna (Pico Veleta, Spain) 
during February 23-24 1995. Observations and d ata reduc- 
tions w ere described sev era l times elsewhere ( ( Andreani 
1994| ), ( [Klein et al. 199e| ), ( [Cimatti et al. 1998| )) and are 



only briefly outlined here. 

For very steep inverted spectra like these, the central 
effective wavelength of the system is 1.25 mm. The beam 
size is 11 arcsec FWHM and the chop throw was set at 60". 

The average sensitivity for each channel, limited mainly 
by atmospheric noise, was of 60 inJy/\/Hz. The sky noise 
was reduced by a large factor using the correlation among 
signals from the channels surrounding the central beam. 

Each source was observed with integration times in the 
range 4000+12000 s with the standard three beam (beam- 
switching + nodding) technique. Atmospheric transmission 
was monitored by frequent sky-dips, which showed zenith 
opacities of typically 0.2-0.4. Calibration was performed us- 
ing Uranus as primary calibrator and quasars of known flux 
as secondary calibrators. The pointing accuracy was checked 
each hour and the average error turned out to be 3". 

Data were reduced assuming that the target sources are 



point-like, or at least that their extent at m,m, wavelengths 
is smallei Lliaii the size uf Lhe ceuLial beam. The uLliei six 
channels were then exploited to reduced the sky noise in the 
central one, by subtracting from the signal in the central 
channel the average value of the atmospheric level computed 
over the outer six ones, for any given elementary integration. 



2.4 Radio data 

Radio VLA data were used to divide the two quasar pop- 
ulations on the basis of the ratio R between the radio flux 
density at 5 GHz, Ssghz, and the optical flux at 4400 A , 
Fb, the radio-quiet objects having R = ^S"' < 10 (e.g. 
( Kellerman et al. 1989 ), ( Falcke et al. 1996 )). Data are t aken 



from 
from 



( Schneider et al. (1991) 1) . (iSchneidCT et al. (1992)|) and 
(Kukula et al. (1998)D, ([Veron fc Veron (1998)1). 



3 EMISSION MECHANISMS IN 
RADIO-QUIET AGN 

Though there has been an active discussion during the past 
several years about the origin of the AGN continuum at in- 
frared wavelengths (the two alternative interpretations being 
either a synchrotron radiation from the AGN itself, or emis- 
sion by circum- nuclear dust), the more recent tendency was 
to clearly favour the latter, at least for radio-quiet objects. 
Crucial observations in this sense were the comparison of 



the millimetric fluxes with IRAS survey data, showing sub- 
millimeter slopes of the SED with spectral indices in excess 
of 3. 

In principle, synchrotron self-absorbed (SSA) sources 
can have a spectral- index steeper than 2.5 if the electron en- 
ergy distribution is dominated by a steep low-energy electron 
population (de Kool & Begelman 1989; Schlickheiser et al. 
1990). In practise, however, with a non-thermal turnover in 
the FIR (at 100-200^tm), such steep indices are not reached 
until mm-wavelengths and hence the available sub-mm data 
clearly indicate thermal processes as the origin of the IR flux 
in AGN. 

Indeed, dust organized in a toroidal structure surround- 
ing the central engine, as envisaged by unifying models of the 
AGN activity, and with a wide temperature distribution up 
to ~ 1500 K (where grains begin to sublimate), re- radiates 
thermally at A > l/im photons absorbed from the p rimary 
source (Neugebauer et al. 1987, ( Sanders et al. 198£| )). This 
interpretation in terms of dust emission naturally accounts 
for the basic observed features of the long-wavelength spec- 
trum of AGN. In particular, it explains the very steep spec- 
tral indices for AGN in the sub-mm as grey-body emission of 
thermal dust with absorption coefBcients depending on fre- 
quency as a{L') ex i/^~^. It also nicely explains the minimum 
in the SED often observed in AGN at A ~ 1 /im and due, 
in this framework, to the cross-over of the primary or stel- 
lar continuum with the emissi on spectrum from d ust at the 
sublimation temperature (e.g. ( panders et al. 1989 ), Granato 
& Danese 1994). Thermal emission at these sub-mm wave- 
lengths is also supported by the detections of CO emission 
f rom low and hi gh -z radio-quiet quasars (see the rev iew by 
(Barvainis 1997); (^chinnerer, Eckart, Tacconi 1998| 



A more com plex situation is considered by (IRowan 



Robinson (1995) ), who interprets the infrared spectrum of 



radio-quiet quasars (RRQs) as contributed both by an AGN 
component with peak emission at A = 3 — 30^m and a star- 
burst emission at A > 30^m. 

At wavelengths longer than 100 fim in the object's rest- 
frame, it is possible that thermal emission from t he host 
galaxy dom in ates the FIR/mm SED of the quasar ((Danese 



et al. 1992), (Fadda et al. 1998)). Testing this is impossible 



with the poor spatial resolution of current instrumentation 
(which would rather require the resolving power of large 
millimeter arrays in plan). 

3.1 A model of dust emission in quasars 

We then assume in the following that the observed FIR/mm 
emission is due to dust surrounding the AGN. We further 
make the rather coarse assumption that the optical proper- 
ties of this dust are roughly similar for all sources, with a 
standard Galactic composition. 

The dust illuminated by the central nuclear source 
reaches an equilibrium temperature which is a function 
of the intensity of the radiation field (hence of the dis- 
tance from the central source) and of the dimensions and 
chemical composition of the assumed grain compon ents. We 
model the SEDs of radio- weak quasars following (Granato 



and Danese (1994' ). The observational data are fitted with 
model spectra produced by a numerical code which solves 
the radiative transfer equation in a circumnuclear dust dis- 
tribution. 
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The latter step is required since in the torii predicted by 
unified models the dust emission is self-absorbed even in the 
near- and mid-IR. The code is quite flexible in dealing with 
different geometries and composition of the dusty medium, 
the only restriction being axial symmetry. It thus allows a 
wide exploration of the parameter space. 

The inner radius ri„ of the dust distribution is set by 
the grain sublimation condition, at Ts — 1500 for graphite 
and Ts — 1000 for silicates. This translates into the condi- 
tion rin ~ 0.5Vi46 pc, where L46 is the luminosity of the 
primary optical-UV emission in units of W^^ erg s~^. 

The details of the model, as well as the effect of the vari- 



ous fre e parameters , have been widely discussed by (Granatc 
and Da iese (1994) ) and ( Granato et al (1997) ). Here we fo- 
cus on the simplest geometry which is in reasonable agree- 
ment with the available observations, the "flared disc", in 
which the scale height of the torus along the z-axis increases 
linearly with the radial distance 



p{r,Q) — Cr ea;p(— 7I cosO|) 



(1) 



where Oj, < < n—Oh and Qh > is the half opening angle 
of the dust-free cones and 7 = for a torus- like distribution. 
The dust distribution within the disc is taken homogeneous 
{(3 = 0). Model parameters are then (i) the outer radius Vmax 
of the dust distribution; (ii) the absorption Ay along typical 
obscured directions, or the corresponding optical depth r at 
0.3 fim, where Ay — 0.61r; (in) the angle O between the 
torus axis and the line-of-sight or equivalently the covering 
factor / = cos(0h) of the torus; (iv) the value of 7 in eq.(l). 
For most of the sources free parameters in the flt are only 
rmax and Ay, while 7 was set to and the covering factor 
has been fixed to / = 0.7, a value not inconsistent with 
the requirements of unified models of AGN (Granato et al. 
1997) (see table 3). Only two objects (0759+651, 1334-1-243) 
were fitted with a dust distribution completely covering the 
nucleus, thicker at the equator than at the poles (in this case 
7>0 and e^ = 0). 

To the first order, the outer radius r,nax mostly deter- 
mines the broadness in wavelength of the IR bump arising 
from dust re-processing, whilst the second parameter. Ay, 
controls the near-IR slope of the SEDs as observed from 
obscured directions, as well as its anisotropy. Therefore Ay 
is in principle testable from measurements of the SEDs of 
obscured AGN or from the anisotropy of mid-IR emission. 
The covering factor is obviously related to the relative power 
reprocessed by the dust in the IR bump and the primary op- 
tical continuum fiux. 

In the simpler geometry with 7 = the total dust mass 
is given by the product between the assumed disc volume 
V = IttcosO^ and the (assumed constant) dust density 
where k^ is the absorption coefficient per 



unit dust mass and r^ the equatorial optical thickness. This 
gives: 






0.2 (- 



-) COS0h Ay 



(2) 



as the total mass in dust as a function of the basic model 
parameters. 



4 RESULTS 



4.1 Masses and sizes of the dust distributions 

Figure 1 shows the spectral energy distribution as vLi, ver- 
sus v in the object's rest-frame of all sources in the sam- 
ple with either FIR or sub-mm detections. The spectra 
are normalized at the 100 fim datum. Some information 
on the source and best-fitting parameters is reported in 
Table 3. Significant constraints on dust masses and sizes 
are estimated for sources with only upper limits to the 
long- wavelength fiuxes. The thick solid curves in Fig.l are 
the best-fit SEDs, while the dotted line is the broad-band 
optical- mm spectrum of a typical local spiral galaxy, used for 
comparison. We adopted to this purpose the SED of M51, 
with its proper normalization corresponding to an assumed 
distance of 9.6 Mpc (see Silva et al. 1998). 

Clearly, while for local objects the hosting galaxy can 
contribute to some extent at the longest wavelengths, for 
high-z quasars the observed spectra are far in excess of what 
would be expected from a redshifted SED of a normal galaxy. 
This refiects the much increased power of the illuminating 
source with respect to what happens in local spiral galaxies. 

The overall quality of the fits is remarkably good, if we 
consider that essentially two parameters {vmax and Ay, the 
inclination angle not affecting significantly the spectrum 
for most of the sources) are free to model the spectral shape. 

Figure 2 plots the mass Md of the circum-nuclear dust 
distribution derived from our spectral best-fits, as a func- 
tion of the B-band luminosity, Ls. The latter measures the 
energy output of the primary nuclear component. Although 
there is a tendency for the lowest- luminosity objects to have 
lower values of Md, the correlation is generally very poor 
and the scatter very large, in particular there is no depen- 
dence on Lb for Lb > 3 10^^ Lq . This is all but reminiscent 
of a simple scaling of the mass in dust with the quasar lu- 
minosity. 

Figure 3 shows the distribution of the dust mass Md as a 
function of redshift. The down-pointing arrows correspond 
to objects with only upper limits to the long wavelength 
fiuxes. The solid horizontal line marks the typical dust mass 
of a nearby spiral (Andreani & Franceschini 1996). The lower 
panel shows the dust masses computed in a go = 0.5, while 
the upper panel reports the corresponding estimates in a 
low density universe (go = 0.1). Low-z objects present a 
large spread in the values of Md, from 10* to 10* Mq. The 
majority of high mass objects are at redshift larger than 1. 
An obvious bias selects at high redshift only luminous and 
massive objects and the lower right part of the plot in figure 
3 is uncovered. 

Figure 4 plots the best-fit outer radius Vmax of the dust 
distribution versus redshift. Md and Vmax are very tightly 
related in our model, since most of the mass comes from 
the outer, less illuminated, lower temperature regions of the 
dust distribution. Both of these parameters are mainly con- 
strained by the long wavelength {mm or sub-mm) flux data. 
Then the plots in Figs. 3 and 4 display quite a similar be- 
haviour, as expected. 

In both cases the low-z objects show widely dispersed 
values in the parameters {Md ~ 10* — 10* Mq) and rmax 
from 200 pc to less than 3 kpc. For the high-z subsample, 
typical values for Md and Vmax are quite signiflcantly larger, 
with Md in excess of lO'^ and r^ax in excess of 1 kpc. 

Note that the procedure of the dust mass estimation 
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is quite robust; we believe that uncertainties in the model 
parameters affect its value by less than a factor two, and our 
basic conclusions are rather secure. 



4.2 A far-IR to optical color-magnitude diagram 
for quasars 

We report in Figure 5 the ratio of the far-IR (I/f/a) to 
the B-band (Lopt) luminosities as a function of Lboi (in so- 
lar units). Lfir is obtained by integrating the rest-frame 
SED from 1/im to 1 mm, while L^pt is found by integrating 
the spectrum between 1200 and 9000 A . It should be first 
noticed in Fig. 5 the wide range covered by the bolomet- 
ric luminosity, which reflects the enormous range of nuclear 
power. 

Figure 5 shows that there is a rough cquipartition of 
the light directly escaping from the nuclear source - and 
measured by the optical flux - and that absorbed by dust 
and reprocessed in the infrared. A remarkable effect appar- 
ent in the flgure is that, while showing a large spread in the 
bolometric luminosity, our sample objects cover a very nar- 
row interval in the infrared to optical luminosity ratio -jf-^-^ 
between ~ —0.1 and ~ 0.2. Only two objects do not fol- 
low this rule: 0759-1-651 presents a distinct behaviour with 
log{Lhoi) = 13.4 and log{ ^/'"^ ) ~ 1.06 due to extinction 
of the optical light. Note, however, that this object is the 
only one of the sample which was not optically selected but 
it was discovered by IRAS. 0844-1-349 on the contrary (at 



log{Li,oi) = 12.6 and log{ 



L) ~ —0.36) presents a high 



optical luminosity with respect to the FIR one (see the cor- 
responding spectrum in Fig. 1). 

In the unified AGN scheme the color axis, ,^'" , is re- 



lated to the viewing angle of the torus, with more inclined 
objects lying in the upper part of the diagram, because of 
the larger FIR emission. The almost complete lack of ob- 
jects in this region of the diagram emphasizes that optical 
surveys, like the present one, provide a strongly biased cen- 
sus of the quasar population: they almost entirely miss ex- 
tinguished objects. Our model of circumnuclear dust (whose 
correctness is proven by the good match with the observed 
spectra) predicts that an unbiased sampling of the quasar 
population (e.g. a selection at long wavelength where r ~ 0) 
would found roughly 25% of the objects to be highly extin- 
guished quasars (i.e. with Oj, > 70°) . 

The origin of the small scatter in the far-infrared to 
optical-UV luminosity ratio of the optically selected quasars 
is not clear and has to be further investigated. In particular, 
there is no evidence for a population of quasars character- 
ized by very low far-IR emission, i.e., in the present inter- 
pretation, by negligible amounts of circum- nuclear dust. In- 
deed, the distribution of upper limits in Fig. 5 may originate 
from the same parent population of the long-wavelength de- 
tections. The lack of objects with low IR/optical luminos- 
ity ratio may be simply due to the sensitivity limits of the 
present millimeter observations but it may also be indicative 
of a narrow intrinsic distribution in the FIR to optical lu- 
minosity ratios. Extensive surveys with sub-mm telescopes 
(SCUBA, CSO) are needed to solve this ambiguity. 



4.3 Radio-loudness versus radio-quiescence 

Figure 6 compares the radio power to the quasar optical 
luminosity. Filled squares correspond to the radio-loud sub- 
population, the other symbols to the radio-quiet one. The 
radio-loud objects are selected according to their ratio R = 
%^ > 10. 

The radio-quiet quasars show a fairly well defined (not 
far from linear) correlation between the two luminosities, 
in spite of the large fraction of radio upper limits, sugges- 
tive of a common physical mechanism originating the two 
emissions. 

The radio activity shows up as a dramatic increase, by 
several orders of magnitude, of the radio power occurring in 
a fraction of the highest optical luminosity quasars. A very 
similar plot of optical versus radio luminosities, scaled down 
in luminosity, was found by Calvani, Fasano, & Franceschini 
(1989) for radio galaxies. Low- luminosity AGN seem in any 
case to provide an unfavourable environment for the ignition 
of radio activity. 

To understand if radio activity is related to some un- 
usual values of the IR to optical luminosity ratios, we have 
reported in Figure 7 the far-IR to optical versus far-IR to 
radio luminosity ratios. The radio active quasars show the 
lowest , "^""^ values. In fact, radio activity does not entail 

any unusual effects in the —tSSii flux ratio and seems to be 
independent of the optical and FIR properties. 

The origin of the millimetric flux in radio-loud objects 



-'-'radio J-^ opt 

Li 



is finally studied in Figure 8, a plot of 

Only very few radio-loud quasars have ^^'"'" values consis- 
tent with a power-law synchrotron extending from the radio 
centimetric to the mm. In the vast majority of objects the 
mm emission is in excess of the radio one, indicating that 
at wavelengths shorter than 1 mm either a high frequency 
synchrotron self-absorbed, or more likely dust emission is 
overwhelming the non-thermal one. 



5 DISCUSSION 

Evidence in favour of the idea that dust emission processes 
dominate the infrared through millimetric emission spec- 
trum of AGN has steadily increased in the last years. 

We are able to confirm this indication with the present 
optically selected sample, but mostly limited to low-redshift 
quasars with far-IR IRAS detections and upper limits to 
the millimetric fiuxes. In more than 70% of the sources the 
sensitivity of the mm observations combined with the IRAS 
detections (see spectra in Figure 1) are enough to exclude a 
synchrotron component based on the steepness of the con- 
tinuum. 

Unfortunately, at redshift larger than 1, with the excep- 
tion of three quasars, objects with mm or sub-mm detections 
(favoured by the strong K-correction) were not usually de- 
tected by IRAS and a similar comparison with mm and sub- 
mm fluxes was not typically possible. The three exceptions 
are QSO0838-h3555 (at z=1.77), QSO1634-h706 (at z=1.34) 
and the weU studied QS01413-I-117 (at z=2.546). In these 
cases, with the addition of BR1202-725 (at z=4.69) detected 
at 150 and 200 fim (rest frame), a dust emission spectrum 
with its very steep convergence in the RJ low-frequency limit 
provides an extremely good fit the the observed SED. 
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A dditional support for this interpretation may also be 



Zghs 



found i ' we consider the ubiquitous presence of a dip in the 



3.39 in the line of sight of QSOOOOO-26 ((Lu et al 



199e), (Molaro et al. 1996)) 



SEDs of Figure 1 around A = 1 ^m for all quasars with good 
near-IR data (naturally explained by the dust sublimation 
hypothesis), and from the generally good agreement pro- 
vided by the dust re-radiation model with the broad-band 
data. 

With this working hypothesis in mind, we have ex- 
ploited long-wavelength observations of quasars to investi- 
gate pr operties of the dusty interstellar media in local and 



up to v jry distant objects. Original motivation for this was 
to achieve a new independent information channel, whose 
results about the metal enrichment in extremely distant ob- 
jects are to be compared with independent analyses of emis- 
sion and absorption metal lines. 

The key to achieve this was to combine long- wavelength 
observations (both IRAS and mm/sub-mm fluxes and upper 
limits) of thermal dust emission with optical-UV data on the 
primary illuminating spectrum. This was essential to con- 
strain the two fundamental parameters of the dust model, 
i.e. the size of the dust distribution (or the average distance 
of dust from the central illuminating source) and the equi- 
librium temperature of dust grains. The obvious remaining 
problem was due to the poor description of the quasar SED 
allowed by the available data, which imply somewhat de- 
generate sets of solutions in terms of these two fundamental 
parameters, of dust grain compositions, and so on. 

We emphasize, however, that exploitation of all avail- 
able data allows quite a robust and unbiased evaluation 
of the dust mass Md, which should be independent of the 
quasar luminosity and redshift. 

Given the relative simplicity of the physics ruling dust 
emission, the hope was to obtain meaningful information 
from a limited dataset. This section is devoted to discuss our 
findings of the previous section, starting with a preliminary 
discussion about the origin of the quasar's long wavelength 
flux. 



5.1 Intervening galaxies: quasars with DLy^ 

systems 

Let us check first the possibility that the observed mm emis- 
sion in the highest redshift quasars might be due to galaxies 
intervening along the line of sight. Recent deep surveys with 
the SCUBA bolometer array on JCMT (Hughes et al. 1998a; 
Small et al. 1997; Barger et al. 1998) have revealed that the 
millimetric sky is populated at faint (~ mjy) flux limits, 
such that millimetric surveys are confusion limited already 
at a few mJy at 850 nm. This indicates that there is a non- 
negligible chance of mm detection in a random direction at 
these fluxes. A way to check for the presence of gas-rich 
galaxies along the quasar sight-line, producing a millimetric 
signal unrelated to the quasar, is offered by absorption line 
studies in the quasar continuum, in particular looking for 
damped Ly^ systems which are believed to correspond to 
already assembled galaxies. 

High resolution spectra available for some objects in 
our sample show indeed the presence of such line systems. 
We consider here the following four: Zats = 4.3829 in the 
line of sight of BRI 1202-0725, Zats = 1.7764 in the line of 
sight of QSO1331-H704, Zats = 4.08 in QSO2237-0607, and 



We have used the metal abundances of the absorber es- 
timated from the absorption lines as an estimate of the gas 
metallicity. The amount of the corresponding dust mass is 
then inferred and compared with the one required to fit the 
FIR/sub-mm emission. This implies the assumption that the 
metallicity of a system is related to its dust content, a rela- 
tionship observed in the Galaxy and i n nearby objects (e.g 
Issa et al. 1990; Sodroski et al., 1995; (Andreani & Frances- 



chini 19961 )). For high redshift objects. Pel et al. (1991) com- 



puted the ratio k between the B-band optical depth, tb, 
and the hydrogen column density, Nh, (with respect to the 
Galactic value kc), which gives an estimate of the dust-to- 
gas mass ratio, as a function of the metallicity. We estimate 
in this way the parameters k/ko, tb {k NH{cm~^)/l(P^) 
and Mdust (see values in Table 4). Mdust can be estimated 
from the k value in the framework of a dust model. Using 
the grain model mentioned in §3.1 we can infer the mass in 
g per H-atom, SRs^s, from: 



4 3 

Jtsys ^ ■Jimod /^obs/ ^mod ^ ^^gr/Hatoiji ^^^ Pgr 

where Ng,./Hatom is the number of grain per H-atom. The 
estimated dust mass becomes: 



Md=NH A ^sys 

where A is the projected area of the Galaxy. The dust 
masses, listed in column 8 of Table 4, turn out to be very 
small even if the dust is spread over a diameter of 10 kpc 
and would produce a negligible mm emission. 

These results are in agreement with the analysis by 
Franceschini and Gratton (1997). Absorbing systems in the 
quasar line of sight with escape velocities larger than 5000 
km/s are not associated with the quasar, and have typi- 
cal metallicities much less than solar (in the range Z/Z© ~ 
0.1-^0.001), in marked contrast with the associated systems 
(v < 5000 km/s) which have metallicities Z/Z© ~ 1 or 
higher. Note that only for a solar metallicity {k/kc ~ 1) the 
amount of dust mass in the absorber would be of the order 
of IO^Mq, close to that in present-day spiral disks, but still 
significantly less than that inferred from the far-IR fiux. 

Altogether, the observed mm emission is very likely due 
to dust associated with the quasar itself. 



5.2 Dust in the quasar and in the host galaxy 

Some limitations of the present analysis have to be consid- 
ered before drawing conclusions. 

The first is that the observational SEDs are interpreted 
assuming that most of the IR- mm fiux comes from dust il- 
luminated by the AGN itself and that the dust in the fiared 
disc is homogeneously distributed from the inner sublima- 
tion surface up to kpc-scale distances. 

For most of the local objects only upper limits to the 
mm fluxes are available. So, it cannot be excluded that for 
some of these we are missing some contributions of cold dust 
from the hosting galaxy at A > 100/xm (most of the mm ob- 
servations were collected with large antennas fed with single 
channel bolometers, having therefore a field of view limited 
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to a fe w arc-sec) . Indeed, Danese et al. (1992) and (Rowan- 



Robinson (1995)) already suggested that the FIR spectrum 



at A > 60fim in low-luminosity AGN may be contributed 
by cold dust in the surrounding galaxy located at large ra- 
dial distances from the AGN and with dust masses larger 
than 1-^3 10^ M©. The relative contribution at far-infrared 
wavelengths of the nuclear and extended component, and 



chini & Gratton, 1997) , rather than triggering and preceding 
it (Silk & Rees, 1998). 

The time evolution seen in Figs. 3 and 4 does not seem 
easily interpretable in terms of a unique class of objects. 
Where all of the observed dust mass at high z has gone? 
Was it swept away from the quasar? Was it consumed to 



form stars in a late phase of star formation? ((Eales & Ed- 



hence ilso the total dust mass, could only be settled by munds 1996)). Note that the difference between low and 



high resolution mapping with planned large mm arrays (e.g. 
LSA/MMA). 

The SEDs of high-z quasars are often poorly con- 
strained, due to the fact that IRAS was not sensitive enough 
to detect them (a few exceptions appear in Fig. 1). In any 
case, comparison of our model with the data suggests the 
existence of a large gaseous structure around the nucleus 
(with 1-5 kpc radius), including a dust mass of up to 10*' 
M0, and a corresponding gas masses of typically Mg = 10^° 
to 10" Mq. 

The mass of the dust distribution around quasars shows 
a very weak dependence on the quasar optical luminosity 
Lb (Fig. 2), and consequently on redshift (Fig. 3). Quasars 
with the lowest Lb tend to have less circumnuclear dust, but 
with a large scatter. Indeed, the dependence of the circum- 
nuclear dust mass on redshift does not seem to follow the 
time evolution of the quasar luminosity function (though it 
is significantly affected by the Malmquist bias). This is also 
illustrated by Fig. 2, where the dust mass shows a very weak, 
certainly not linear, dependence on nuclear luminosity. 

Concerning the overall extent rmax of the dust distri- 
bution. Fig. 4 reveals that local quasars and Seyfert galaxies 
(at z < 0.2) show dust structures out to a maximum radius 
of 1 kpc. For higher redshift quasars Vmax increases up to 
values of several kpc, indicating that these structures are 
not merely circum-nuclear dust torii, as seen in local Seyfert 
galaxies, but have rather scale-lengths comparable or even 
larger than those of the bulge of the putative host galaxy. 
How reliable are these estimates? We claim that whenever 
the IRAS flux limits, combined with mm detections, mean- 
ingfully constrain the model SED, there is no margin to 
decrease rmax by increasing the dust temperature Td and 
reducing the mass Md, because this would imply a strong 
and unacceptable violation of the IRAS limit at typically 60 
to 25 fim. Certainly in all these cases our estimate of Vmax 
should be quite reliable. 

The result about r,„ax is indicative of a medium heavily 
enriched by dust and metals on a quite large scale around the 
central AGN, as would be expected for an already large and 
massive galaxy in an advanced stage of chemical evolution. 

At the sensitivity limits of the present millimetric obser- 
vations there is no evidence for a sub-population of quasars 
with no or significantly depressed far-infrared emission. As 
illustrated by Fig. 5, the far-infrared emission scales typi- 
cally linearly with the optical quasar luminosity. Dust seems 
to be a ubiquitous presence at any redshifts, whenever a sen- 
sitive enough observation has targeted it. 

Though needing confirmation by further sensitive ob- 
servations with the new-generation sub-millimetric imagers 
(SCUBA, CSO), this result, also in combination with metal 
line observations, tends to provide support to the view that 
the quasar phenomeno n follows the formation of a large 
spheroidal galaxy (e.g. (Hamann fc Ferland 1995), Frances- 



high redshift objects would be further enhanced in an open 
universe: dust masses at redshifts larger than 2 for Q ~ 0.01 
are higher by a factor between 2 and 5. 

The alternative possibility is that we are observing at 
high- and low-redshifts two distinct quasars populations. En- 
vironmental conditions at high redshifts could favour the 
formation of higher mass black holes, while low-z AGN may 
be related to the refueling of small dead black holes in later 
type galaxies, where gas is still available (Haehnelt & Rees 
1993; Cavaliere & Vittorini 1998). 

We finally note that our optically-selected sample has a 
clear bias against optically extinguished quasars with the 
torus viewed edge-on, as we have shown in Sect. 4.2. In 
addition, an entire early phase in which the object is still 
obscured by the products of massive galaxy-forming star- 
bursts could have escaped detection. Only large-area surveys 
at long-wavelengths, such as those to be performed with the 
ISO, SIRTF and FIRST missions, will be able to discover 
them. 



6 SUMMARY AND CONCLUSIONS 

We have performed a detailed analysis of broad-band spec- 
tra for a rich unbiased sample of optically-selected quasars 
with far-IR and mm data, with an emphasis on the portion 
of the SED dominated by dust thermal emission. The ex- 
ploitation of flux data over a wide frequency interval has 
allowed us to successfully constrain the basic parameters of 
the dust distribution around the quasar, i.e. the size rmax, 
temperature Td and total dust mass Md. 
Our main results are summarized below. 

• We flnd that, when detailed enough data exist, a 
model of emission by circum-nuclear dust illuminated by 
the optical-UV flux of the quasar provides a very good fit to 
the observed SED, including the steep sub-millimetric con- 
vergence and the observed minimum at around A = 1 /im. 
This supports the concept of a thermal origin for the IR 
SED in radio quiet quasars and tends to rule out alternative 
interpretations. 

Radio activity in quasars does not seem to entail unusual 
properties of far-IR and optical emission with respect to the 
radio-quiet population. 

• The mass of dust around high redshift {z > 1) quasars 
is typically observed to fall between lO'' and 10^ M©, with 
no appreciable dependence on z over a very large interval 
from z — 1 to z — 5. Similarly, the radial size of the dust 
distribution rmax runs from rmax — 1 kpc to ~ 10 kpc over 
the same redshift interval. These results indicate the pres- 
ence around quasars of massive dust distributions on scales 
comparable with those of the hosting galaxy. The latter have 
to be already large systems observed in advanced stages of 
chemical evolution during or after a very active phase of star 
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and metal production. This evolutionary pattern looks fun- 
damentally different from that expected for disk-dominated 
galaxies, and requires an early accelerated phase of star- 
formation, probably related to the formation of a massive 
spheroid. 

• At the sensitivity limits of the present data, the ra- 
tios of the bolometric far-infrared to optical-UV luminosi- 
ties show a narrow distribution, implying a roughly con- 
stant fraction of optical-UV light being reprocessed into the 
IR, with no evidence for a population of metal-poor opti- 
cally dominated objects. Although more sensitive mm sur- 
veys (e.g. with SCUBA) are needed to better probe the faint 
end of the distribution, this tends to support the view that 
the quasar phase parallels and/or follows the formation of a 
large hosting galaxy, rather than preceding it. 
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Figure captions 

Figure 1 Spectral energy distributions (S.E.D.s) of the 
quasars of the sample. Each panel reports the name and the 
redshift of the object. The solid line represents the predicted 
S.E.D. according to the model described in the text (§3). The 
dashed-line is the SED of a template nearby dusty spiral 
with a dust mass of 10^ Mq. 

Figure 2 The dust mass computed according to the 
model prescriptions described in the text (§3.1) as a function 
of the B-band luminosity. Down arrows correspond to those 
sources with only upper limits at FIR/sub-mm wavelengths. 

Figure 3 The dust mass as a function of redshift in 
an open universe (go = 0.1, A = 0) (upper panel) and in 
a flat universe (go ~ 0.5, A = 0) (lower panel). Down ar- 
rows correspond to those sources with only upper limits at 
FIR/sub-mm wavelengths. Note how the difference between 
the nearby and the distant source masses increases for an 
open universe. 

Figure 4 The maximum radial size, rmax, of the torus 
surrounding the central source for the quasars of our sample. 



Figure 5 A colour-luminosity log{- 



Ltoi plot for 



the quasars of the sample. Filled circles correspond to low- 
z (z < 1) objects. Open squares to intermediate redshift 
(1 < 2: < 2) sources, asterisks to highz (z > 2) ones. The two 
outliers are 0759-^651 (at LpiR/Lopt = 1.06) and 0844+349 



(at log{- 



-0.36). 



Figure 6 Radio power as a function of the B-band ab- 
solute magnitude. Filled squares correspond to radio-loud 
objects, while open circles to radio-weak sources. Down ar- 
rows show those objects observed at 5 or 1.4 GHz but not 
detected, hence assumeed to belong to the radio- weak pop- 
ulation. 

Figure 7 A colour-colour plot showing the behaviour 
of the rest-frame luminosity ratio, I/ioonm/iopt, versus 
Lioofim/Lradio- Filled syuibols correspond to the radio- loud 
population, while open ones to the radio-weak population. 

Figure 8 A colour-colour plot showing the behaviour 
of the rest-frame luminosity ratio, Limm/Lradto, versus 
Limm/Lopt- Symbols as in figure 7. 
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TABLE 1 
Photometry of radio-quiet quasars 



Name 



m_B 



mv 



m_R 



'H 



'K 



Fi 



2^m 



25jUr7i 



Feo^ur^ 



FlOOi, 



FsSOi, 



Fi^ 



M 



B 



0000-26 


4.11 






17.53 


0019-15 


4.52 






19.00 


0027-F05 


4.21 






21.93 


0046-29 


4.01 






19.30 


0050-M2 


0.06 


14.41 


14.03 


13.63 


0051-27 


4.40 






20.18 


0054-M4 


0.17 




15.71 


15.40 


0101-30 


4.07 






20.06 


0103-FOO 


4.43 






18.70 


0104-F02 


4.17 




19.73 




0114-08 


3.16 




20.50 




0151-00 


4.19 






18.65 


0157-FOO 


0.16 


16.15 


15.69 




0241-01 


4.04 






18.40 


0245-06 


4.20 






18.20 


0245-00 


2.11 


18.68 


18.41 


17.90 


0247-00 


1.45 


17.95 


17.87 


17.70 


0247-02 


1.93 


18.70 


18.20 


18.20 


0249-F02 


2.81 


18.86 


18.49 


18.30 


0250-FOl 


2.63 


18.85 


18.62 


17.90 


0252-FOl 


2.45 


18.23 


18.07 


17.20 


0258-F02 


2.52 


18.00 


17.85 


17.70 


0259-FOl 


1.74 


18.62 


19.53 


18.20 


0307-F02 


4.37 






20.39 


0344-F02 


3.37 






20.43 


0345-FOl 


3.63 






19.94 


0351-10 


4.36 






18.70 


0710-F45 


0.05 


14.57 


14.27 


14.11 


0751-F56 


4.28 




21.90 


19.91 


0759-F65 


0.14 




15.50 


14.30 


0838-F35 


1.77 




16.00 


16.10 


0838-F77 


0.13 


16.30 


15.70 




0842-F34 


2.12 




17.00 


17.00 


0844-F34 


0.06 


14.83 


14.50 





9.91 14 



3.69 5.25 10 



72 



30 9 



37 



43 



<150 
<120 
<120 
<150 
52.4 549±11 



45 



23 24 



2.29 3.09 4.79 
7.94 9.55 12.58 



24 



26 



<120 

72±7 

< 90 

< 75 
<120 
<120 
<150 

137±38 

< 75 

<120 

< 90 

< 90 

< 60 

< 75 

< 90 

< 60 
<105 

< 75 

< 60 

< 60 

< 90 
233±27 

< 70 
360±40 
100±30 

34±5 

< 90 
26±29 



<140 
<160 
<150 
<135 
1097±20 

<100 
86±13 
<100 
<150 
<210 
<120 
<180 
520±58 
150±40 

<105 

< 90 

< 90 
<120 

< 90 
<105 
<120 
<120 
<150 

< 90 
<120 

< 60 
551±33 

<100 
610±50 
150±50 
103±7 
<200 
204±40 



<120 
<150 
<120 
<120 
2293±17 

<120 
291±9 
<135 
<150 
<160 
<135 
<150 
2377±56 
<120 

<120 
<100 

< 90 
<105 

< 75 
<135 
<150 
<240 
<120 
<180 
<100 

< 60 
864±44 

<180 

1750±60 

200±40 

174±9 

<120 
163±41 



<375 
<480 
<370 
<310 

2959±51 

<270 
794±26 
<330 
<500 
<360 
<400 
<360 
2322±130 
<300 

<360 

<420 

<420 

<320 

<345 

<700 

<300 

<700 

<700 

<700 

<700 

<360 
1439±107 

<600 

1810±114 

800±100 

426±30 

<390 
294±97 



<27.0 



18.4±4.4 
(225±75)* 



20.9±5.2 



<27.0 
<16. 
<33.0 
<19. 



<8.0 

<14.0 
<4.3 
<3.6 

<4.2 
<3.0 
<3.6 
<3.3 
<4.5 
<6.0 
<4.0 
<2.7 
<6.5 
<3.6 
<13.0 
<6.2 
<6.2 
<12.7 
5.9±2.7 
<6.5 
<8.8 
<14.3 
<3.0 
<4.0 
<5.4 

<2.1 

<4.0 

42* 

4.5±1.5 

< 3.0 
4.1±1.3 

<3.4 



<0.3" 

<0.3 
3.1±0.2 

< 0.3 

(1.8 ±0.4)" 
<0.2 

<0.3 



7.9±0.3 



<0.13 
<0.15 
<0.12 
<0.08 

<0.1 
<0.10 

<0.3 
<0.08 
< 0.2 

<0.3 

<0.2 

<0.3" 

1.8±0.1 

<0.6 
12±0.5 

<0.4 



-28.80 
-27.60 
-23.96 
-26.96 
-23.36 

-27.16 
-24.11 
-26.27 
-27.60 
-26.61 
-28.68 
-27.40 
-23.81 
-27.80 
-27.50 
-27.01 
-26.93 
-27.17 
-27.40 
-27.40 
-27.77 
-28.03 
-26.64 
-26.44 
-25.86 
-26.58 
-27.80 
-22.46 
-26.70 
-24.23 
-28.71 
-23.73 
-28.61 
-23.07 



TABLE 1 — Continued 



Name 



m_B 



mv 



m_R 



'H 



'K 



Fr2i_ir. 



25jUr7i 



Feoi, 



FioOl, 



FsSO^ur^ 



Fi^ 



M 



B 



0856+17 


2.32 






19.00 


18.70 


0903+17 


2.77 






18.00 




0906+48 


0.11 


16 


46 


16.10 




0910+56 


4.03 






22.10 


20.87 


0914-62 


0.05 


14 


07 


13.55 




0951-04 


4.35 








19.00 


0952-01 


4.43 








18.70 


0953+47 


4.45 






21.37 


19.47 


0956+12 


3.30 






17.60 


17.40 


1001+05 


0.16 


16 


38 


16.23 


16.20 


1013+00 


4.38 








19.10 


1017+10 


3.15 






18.40 


17.30 


1033-03 


4.50 








18.50 


1114-08 


4.50 








19.70 


1117-13 


3.96 








18.10 


1126-04 


0.06 


14 


97 


14.61 


13.50 


1144-07 


4.14 








18.80 


1158+46 


4.73 








20.21 


1202-07 


4.69 








17.50 


1208+10 


3.80 








17.50 


1211+14 


0.08 


14 


40 


14.19 


13.90 


1212+08 


2.35 






18.00 


18.20 


1212+14 


1.64 






17.90 


17.10 


1230+16 


2.70 






17.40 


17.60 


1231+29 


2.01 






16.00 


16.20 


1235+08 


2.90 






17.80 




1246-05 


2.21 


17 


06 


16.70 


16.20 


1247+26 


2.04 






15.60 


15.20 


1247+34 


4.89 






22.60 


20.40 


1301+47 


4.00 






22.20 


21.31 


1302-14 


4.04 








18.40 


1308-02 


2.85 






18.70 


18.40 


1308-01 


2.55 






18.10 


17.60 


1328-04 


4.20 








19.00 



10 



34 



47 14 



16 



79 26 



75 



30 



< 90 


<150 


<100 


<360 


< 70 


<120 


<105 


<360 


39+6 


87+8 


172+10 


291+32 


<75 


< 60 


<105 


<300 


490+30 


1170+99 


2400+99 


1820+120 


< 75 


<210 


<110 


<450 


110+30 


120+50 


1100+40 


2470+160 


< 90 


<100 


<110 


<600 


<120 


<120 


< 90 


<600 


31+6 


36+12 


27+9 


<69 


<150 


<160 


<150 


<420 


< 75 


<200 


< 60 


<270 



< 90 
<150 

<180 

172+38 

<150 

< 75 

< 90 
<120 
<180 

< 90 
<170 

< 90 

< 90 

<240 

< 90 
<120 



< 60 
<150 

<150 
362+51 
<150 
<150 
<150 
<120 
<150 
<240 
<120 

< 90 

< 60 

<150 
<150 
<300 



<150 
<220 

<90 
305+53 
<100 
<120 
<150 
<120 
<120 
<150 
<105 
<100 
< 90 

<150 
<120 
<120 



12+4 



104+19 309+35 669+26 1172+134 



<300 
<230 

<270 
689+119 
<450 
<300 
<360 
<300 
<330 
<300 
<240 
<240 
<450 

<300 
<330 
<300 



50+7 
(92+38)* 

<32 



<4.8 

<6.0 

<2.8 

<6.0 

<6.0 

<3.0 

2.8+0.6 

<2.7 

<6.0 

<3.0 

<7.5 

3.7+1.2 

3.5+0.7 

4.1+0.8 

4.1+0.8 

<3.2 
5.9+1.0 

<6.0 
12.6+2.3 

4.1+0.8 
<2.8 
<6.0 
<9.6 

7.5+1.4 
<6.3 
<4.0 
<6.0 
<6.0 
<1.5 
<5.3 
<3.5 
<6.0 
<7.0 
<3.3 



1.92+0.18 

<0.3 
<0.3 
16+1. 



<0.3 
<0.2 

<0.3" 

0.5" 

<0.3 
<0.6" 

0.8" 



1 

0.72+0.08 

<1 



-26.38 
-27.92 
-22.74 
-25.40 
-23.60 
-27.20 
-27.70 
-27.30 
-29.13 
-23.51 
-27.40 
-27.98 
-27.57 
-26.70 
-28.10 
-22.77 
-27.50 
-26.90 
-28.62 

-28.10 
-24.04 
-27.37 
-26.66 
-28.44 
-29.01 
-28.27 
-28.02 
-29.39 
-27.30 
-24.90 
-27.80 
-27.39 
-28.50 
-27.20 



TABLE 1 — Continued 



Name 


z 


ms 


mv 


m_R 


Fj 


Fh 


Fk 


Fl2^^.m 


F25^^.m 


Feodum 


FlOO^um 


FsSO^um 


Fl,25mm 


Fecm 


Mb 


1331-01 


1.88 




17.90 


16.70 


















<7.0 


0.9±0.1 


-26.86 


1334-F24 


0.10 




15.10 


14.30 








680±70 


<820±80 


770±70 


550±50 




4.4±1.1 


<6 


-24.04 


1335-04 


4.40 






19.10 


















10.3±1.0 




-27.30 


1346-03 


4.01 






19.40 


















<4.2 




-26.80 


1351-F64 


0.08 


14.54 


14.28 


13.70 


6.61 


7.08 


10.9 6 


173±5 


600±6 


810±8 


1350±26 




<2.5 


13.0±0.4 


-24.04 


1402-F43 


0.32 




16.50 


14.70 








190±10 


290±10 


590±50 


1150±100 


106* 






-24.82 


1411-F44 


0.08 




14.99 


13.90 


5.62 


8.32 


17.38 


115±10 


160±12 


162±17 


<175 


<36 


3.9±1.0 


<0.6" 


-23.64 


1413-Ml 


2.54 




17.00 


16.20 








< 93 


<126 


230±38 


370±78 


44±8 
(189±56)'' 


18±2 
(224±30)* 


(1.95 ±0.13) 


-28.60 


1426-01 


3.42 


18.44 


17.80 


17.50 


















<6.6 




-28.40 


1426-FOl 


0.08 


15.21 


14.87 


13.60 


5.89 


8.71 


15.85 


124±33 


171±46 


318±47 


<315 




<4.3 


0.84±0.08 


-23.34 


1435-F38 


1.60 




17.00 


17.30 


















<12.0 


70 


-27.48 


1440-F35 


0.07 


14.90 


14.58 


14.28 


10.0 


14.13 


25.12 


98±13 


208±15 


652±21 


1061±64 




<3.5 


0.8±0.1 


-23.35 


1449-F58 


0.21 




16.00 


15.30 








43±9 


60±11 


123±16 


34252 




<3.1 




-24.38 


1500-F08 


3.96 






18.84 


















<4.0 




-27.10 


1543-F48 


0.40 


16.05 


16.65 


15.70 








58±15 


126±18 


348±26 


485±79 




<4.1 


1.04±0.09 


-25.08 


1548-F46 


3.54 




19.62 


19.24 


















<1.5 


<0.3 


-26.96 


1613-F65 


0.12 


16.03 


15.49 




4.57 


6.31 


10.72 


87±12 


231±14 


635±19 


1090±59 




<3.1 


3.0±0.2 


-23.43 


1634-F70 


1.33 


14.90 


14.66 


14.60 


6.61 


9.55 


6.76 


61±13 


147±14 


318±23 


343±72 


<47 


<4.2 


l.liO.l 


-29.60 


1640-F46 


3.70 




20.06 


19.52 


















<3.4 


<0.3 


-26.80 


1643-F46 


3.79 




20.87 


20.26 


















<3.0 


<0.3 


-26.10 


1643-F46 


3.83 




21.29 


20.61 


















<2.4 


<0.3 


-25.80 


1644-F47 


3.69 






20.69 


















< 2.9 




-25.01 


1700-F51 


0.29 


15.43 


15.10 


14.70 


4.37 


6.46 


12.6 


110±23 


270±21 


420±36 


582±88 




<4.4 


7.2±0.2 


-26.04 


1700-F64 


2.72 


16.30 


16.13 


15.70 








< 36 


< 40 


56±19 


120±50 




<6.0 




-29.50 


1821-F64 


0.29 


14.20 


14.20 


13.90 








236±19 


373±18 


953±32 


2164±82 


<45* 


<3.7 




-26.94 


2041-10 


0.03 


14.31 


14.00 


13.46 


14.72 


24.49 


41.01 


350±31 


667±44 


1508 44 


1607 104 




4±1 


5.5±0.2 


-22.76 


2047-FOl 


3.80 




20.54 


19.71 








<105 


<100 


<150 


<400 




1.9±0.5 


<0.3 


-26.90 


2048-FOl 


3.21 




20.54 


20.65 








< 75 


< 90 


<120 


<550 




<1.5 




-24.76 


2112-F05 


0.45 


16.06 


15.77 


15.40 


2.88 


3.80 


6.61 


71±13 


73±25 


105±19 


<177 




<2.1 


0.9±0.1 


-25.95 


2130-F09 


0.06 


14.92 


14.64 


13.50 


7.90 


14.13 


25.12 


188±8 


378±14 


480±13 


<1888 




2.1±0.7 


2.1±0.1 


-22.98 


2132-FOl 


3.19 




19.67 


19.78 








<90 


<120 


<210 


<330 


<12.0 


<2.4 


<0.2 


-26.17 


2203-F29 


4.40 






20.87 








<75 


<120 


<150 


<270 




<4.8 




-25.96 


2212-16 


4.00 






18.60 


















<6.6 




-27.60 


2214-M3 


0.06 


15.29 


14.90 


14.64 


7.9 


13.2 


12.6 


61±7 


95±12 


337±11 


<282 


<49. 


3.5±0.6 


<0.2 


-22.98 



TABLE 1 — Continued 



Name 



ms 



mv 



m_R 



'H 



' K 



-T 12jUr7i -T 25jUr7i -^ 60jUr7i -^ lOOjUm 



- SSOjUm 



Fi 



25mm 



6cm 



Mf 



2235-03 


4.23 


2237-06 


4.55 


2248-12 


4.16 



18.60 
18.30 





3.0 


<0.6" 


-27.80 


<21.0 


<4.0 


<0.6" 


-28.10 


<30.0 




<0.6" 





^radio measurements at 20cm 
'^measurements at 450 fim 
■^measurements at 1.1mm 
■^measurements at 350 /im 

Note. — All fluxes are given in mjy 



TABLE 2 

Photometry of radio-loud quasars 



Name 



m_B 



mv 



m_R 



'H 



'K 



Fi 



2jj,m 



2hjj,m 



^QOjir. 



FlOOi. 



850jum 



Fi^ 



M 



B 



0251-00 


1.67 


18.59 


18.49 


16.90 


0256-00 


3.36 


18.72 


17.50 


18.10 


0256-00 


1.74 


18.76 


18.68 


18.00 


0923-F39 


0.69 


17.92 


17.86 


17.10 


1050-00 


4.29 






18.80 


1331-M7 


2.08 


16.84 


16.71 


16.10 


1343-F38 


1.84 




18.50 


18.00 


1348-F39 


1.59 






19.00 


1353-M8 


0.05 


15.01 


14.22 




1417-F38 


1.83 




19.30 


18.50 


1500-F04 


3.67 






18.00 


1508-F57 


4.30 






18.90 


1557-F03 


3.90 






19.80 


1745-F62 


3.89 






18.80 


2331-F02 


4.09 




21.13 


19.98 






<120 


< 70 


<120 


<300 


77 


145±27 


<120 


150±42 


<400 




<120 


<240 


<180 


<700 


22 











41 



< 90 



<150 



<150 



< 50 
<120 



< 40 
<180 



< 90 
<165 



<300 



482±35 1380±37 2218±48 1986±125 



<25.0 
<30.0 



<150 
<450 



<8.5 

<11.8 

138±9 

3744±53 

<9.0 

<4.5 

<5.4 

39.7±9.3 

3.1±0.7 

65.6±3.6 

<7.2 

(38.0±6.0)<^ 

23.0±4.5 

50±7 

<3.6 



12.9±0.1 
3.73±0.04 
506.6±0.2 

7570 

; 10.6±0.2)'' 

713±97 

416 

89 

100±6 

871 
169±17 
(282±29)'' 
690±19 
580±55 
3.3± 0.1 



-26.57 
-28.58 
-27.96 
-24.92 
-26.90 
-28.20 
-26.28 
-25.49 
-22.37 
-25.49 
-27.62 
-27.03 
-25.95 
-26.95 
-26.41 



^radio measurements at 20cm 

■^measurements at 1.1mm 

Note. — All fluxes are given in mjy 



TABLE 3 

Dust masses and sizes 



name 


^miTi 


(pc) 


r 


e. 


Mdu.t 

(Me) 

(fi = l) 


Mdu.t 

(Me) 

(fi = 0.2) 


0000 -26 


< 800 


2.48 


30 


0° 


<5.2 10^ 


<1.5 10® 


0019-1522 


< 800 


1.58 


40 


45° 


<5.2 10^ 


<1.6 10® 


0027 -F0521 


< 8000 


0.39 


60 


45° 


<4.8 10^ 


<1.4 10^ 


0046 -293 


< 800 


1.23 


30 


45° 


<2.3 10^ 


<6.5 10^ 


0050 -M24 


2000 


0.24 


30 


45° 


5.9 10^ 


6.0 10^ 


0051 -279 


< 800 


0.92 


40 


45° 


<1.7 10^ 


<5.2 10^ 


0054 -M44 


1500 


0.30 


50 


0° 


8.5 10^ 


9.1 10^ 


0101-304 


< 1000 


0.88 


30 


45° 


<1.9 10^ 


<5.3 10^ 


0103 -F0032 


< 400 


1.71 


30 


45° 


<1.2 10^ 


<3.4 10^ 


0104-F0215 


< 400 


1.71 


30 


45° 


<1.2 10^ 


<3.3 10^ 


0114-0857 


< 2000 


0.65 


30 


45° 


<4.2 10^ 


<1.1 10® 


0151 -0025 


< 800 


1.54 


30 


0° 


<3.7 10^ 


<1.1 10® 


0157 -FOOl 


2000 


0.40 


60 


45° 


3.0 10^ 


3.3 10^ 


0241 -0146 


< 1000 


1.64 


30 


0° 


<6.6 10^ 


<1.8 10® 


0245 -002 


< 2000 


1.12 


30 


45° 


<1.2 10^ 


<2.4 10® 


0247 -002 


< 2000 


1.06 


30 


45° 


<1.1 10® 


<1.8 10® 


0247 -023 


< 1000 


1.11 


50 


45° 


<4.9 10^ 


<9.3 10^ 


0249 -F022 


< 1000 


1.53 


30 


45° 


<5.6 10^ 


<1.3 10® 


0250 +0U 


1000 


1.28 


60 


45° 


8.0 10^ 


1.8 10® 


0251 -000 


< 1500 


0.94 


60 


45° 


<9.8 10^ 


<1.7 10® 


0252 +on 


< 600 


1.51 


30 


45° 


<2.0 10^ 


<4.2 10^ 


0256 -000 


< 600 


2.59 


30 


45° 


<5.7 10^ 


<1.4 10® 


0256 -003 


< 1000 


0.83 


30 


45° 


<1.7 10^ 


<3.0 10^ 


0258 -F021 


< 600 


1.62 


30 


45° 


<2.2 10^ 


<4.8 10^ 


0259 +0U 


< 2000 


0.73 


60 


45° 


<1.2 10® 


<2.2 10® 


0307 -F0222 


< 1000 


2.60 


30 


45° 


<1.7 10® 


<4.9 10® 


0344 -F0222 


< 2000 


0.62 


30 


45° 


<3.8 10^ 


<9.6 10^ 


0345 -F0130 


< 1000 


0.85 


30 


45° 


<1.7 10^ 


<4.6 10^ 


0351 -1034 


< 800 


1.71 


30 


45° 


<4.6 10^ 


<1.4 10® 


0710 -F457 


1000 


0.17 


30 


45° 


6.9 10^ 


7.0 10^ 


0751 -F5623 


< 1000 


0.84 


30 


45° 


<1.7 10^ 


<5.0 10^ 


0759 -F651 


3000 


0.40 


80 


a 


3.6 10^ 


3.9 10^ 


0838 -F3555 


400 


3.15 


30 


45° 


3.9 10^ 


7.0 10^ 


0838 -F770 


2000 


0.19 


60 


0° 


7.1 10^ 


7.5 10^ 


0842 -F3431 


800 


1.97 


60 


0° 


1.2 10® 


2.4 10® 



TABLE 3 — Continued 



name 




'^'min 

(pc) 


r 


e. 


Mdtist 

(Me) 

(fi = l) 


Mdtist 

(Me) 

(fi = 0.2) 


0844 +349 


600 


0.16 


10 


10° 


3.8 10^ 


3.9 10^ 


0856 +1713 


< 1000 


0.91 


60 


0° 


<4.0 10^ 


<8.2 10^ 


0903 +1734 


< 400 


1.68 


60 


0° 


<2.2 10^ 


<4.9 10^ 


0906 +484 


2000 


0.17 


60 


0° 


5.5 10^ 


5.8 10^ 


0910 +5625 


< 4000 


0.53 


60 


0° 


<2.2 10^ 


<6.2 10^ 


0914-6206 


800 


0.26 


60 


0° 


2.0 10^ 


2.1 10^ 


0923 +392 


400 


0.56 


30 


45° 


1.0 10^ 


1.3 10^ 


0951 -0450 


< 800 


1.34 


30 


45° 


<2.8 10^ 


<8.3 10^ 


0952-0115 


400 


1.79 


30 


0° 


1.3 10^ 


3.8 10^ 


0953 +4749 


< 1000 


1.09 


30 


45° 


<2.9 10^ 


<8.7 10^ 


0956 +1217 


< 200 


2.53 


30 


0° 


<6.2 10^ 


<1.6 10^ 


1001 +054 


100 


0.20 


30 


0° 


9.7 10^ 


1.0 10^ 


1013 +0035 


< 800 


1.33 


30 


45° 


<3.0 10^ 


<9.0 10^ 


1017+1055 


1000 


1.49 


30 


0° 


5.5 10^ 


1.3 10® 


1033 -0327 


600 


1.86 


60 


0° 


6.1 10^ 


1.8 10® 


1050 -0000 


< 400 


1.64 


60 


45° 


<2.1 10^ 


<6.1 10^ 


1114-0822 


1500 


1.04 


60 


0° 


1.2 10^ 


3.6 10® 


1117-1329 


600 


1.90 


60 


0° 


6.4 10^ 


1.8 10® 


1126-041 


2000 


0.16 


60 


0° 


5.0 10^ 


5.1 10^ 


1144-0723 


1000 


1.53 


60 


0° 


1.1 10^ 


3.2 10® 


1158 +4635 


< 1000 


0.92 


60 


45° 


<4.1 10^ 


<1.3 10® 


1202 -0725 


800 


2.94 


60 


0° 


2.7 10^ 


8.2 10® 


1208+1011 


800 


2.39 


30 


0° 


8.9 10^ 


2.4 10® 


1211 +143 


2000 


0.23 


30 


0° 


5.2 10^ 


5.4 10^ 


1212 +0854 


< 1000 


1.38 


30 


0° 


<4.7 10^ 


<9.6 10^ 


1212 +1445 


< 2000 


1.08 


60 


0° 


<2.3 10^ 


<3.9 10® 


1230 +1627 


1000 


2.01 


60 


0° 


2.0 10^ 


4.4 10® 


1231 +2924 


< 600 


2.94 


30 


0° 


<7.5 10^ 


<1.4 10® 


1235 +0857 


< 400 


1.94 


30 


45° 


<1.5 10^ 


<3.4 10^ 


1246 -0542 


< 800 


2.31 


30 


0° 


<8.2 10^ 


<1.7 10® 


1247 +2647 


< 800 


3.62 


30 


0° 


<3.2 10^ 


<6.3 10^ 


1247 +3406 


< 400 


1.58 


30 


0° 


<9.7 10^ 


<3.1 10^ 


1301 +4747 


< 2000 


1.88 


30 


45° 


<2.3 10^ 


<6.3 10^ 


1308-0214 


< 800 


1.30 


30 


45° 


<2.6 10^ 


<6.0 10^ 


1308-0104 


< 400 


1.51 


30 


45° 


<3.6 10^ 


<7.7 10^ 



TABLE 3 — Continued 



name 




'^'min 

(pc) 


r 


e. 


Mdtist 

(Me) 

(fi = l) 


Mdtist 

(Me) 

(fi = 0.2) 


1328 -0433 


< 400 


1.50 


30 


45° 


<8.6 10^ 


<2.5 10^ 


1334 +243 


3000 


0.32 


60 


a 


1.3 10^ 


1.4 10^ 


1335 -0417 


2000 


1.33 


60 


0° 


3.8 10^ 


1.1 10^ 


1348 +39 


8000 


0.54 


60 


45° 


1.0 10^ 


1.8 10^ 


1351 -F640 


2000 


0.24 


30 


45° 


5.6 10^ 


5.8 10^ 


1353 -M86 


2000 


0.19 


30 


45° 


3.4 10^ 


3.5 10^ 


1402 -F437 


2000 


0.67 


30 


0° 


4.3 10^ 


4.9 10^ 


1411 +442 


450 


0.17 


17 


0° 


6.1 10^ 


6.1 10^ 


1413 -M17 


400 


3.58 


60 


45° 


1.0 10^ 


2.2 10^ 


1417 -F385 


8000 


1.36 


60 


45° 


1.4 10^ 


2.6 10^ 


1426 -F015 


600 


0.19 


50 


0° 


5.3 10^ 


5.5 10^ 


1440 -F356 


1000 


0.22 


30 


0° 


1.1 10^ 


1.2 10^ 


1449 -F588 


4000 


0.29 


30 


0° 


3.3 10^ 


3.6 10^ 


1500 -F0431 


< 400 


2.08 


60 


45° 


<3.4 10^ 


<9.1 10^ 


1508 -F5714 


< 2000 


1.36 


30 


0° 


<1.8 10^ 


<5.2 10^ 


1543 -F489 


2000 


0.66 


60 


0° 


8.5 10^ 


1.0 10^ 


1557 -F0313 


4000 


1.00 


60 


45° 


7.8 10^ 


2.2 10^ 


1613 -F658 


2000 


0.27 


60 


0° 


1.4 10^ 


1.5 10^ 


1634 -F706 


600 


2.70 


30 


0° 


6.1 10^ 


9.8 10^ 


1700 -F518 


1000 


0.67 


30 


0° 


1.1 10^ 


1.2 10^ 


1700 -F6416 


400 


2.72 


30 


0° 


2.9 10^ 


6.5 10^ 


1745 -F6227 


4000 


1.54 


60 


45° 


1.9 10^ 


5.2 10^ 


1821 -F643 


1000 


1.04 


60 


0° 


5.3 10^ 


6.0 10^ 


2041 -109 


2000 


0.14 


60 


0° 


3.7 10^ 


3.7 10^ 


2047 -F0123 


600 


0.96 


60 


45° 


1.6 10^ 


4.4 10^ 


2048 -F0126 


< 600 


0.59 


60 


45° 


<6.1 10^ 


<1.5 10^ 


2112 -F059 


200 


0.69 


30 


0° 


4.6 10^ 


5.5 10^ 


2130 -F099 


800 


0.17 


30 


0° 


4.7 10^ 


4.8 10^ 


2132 -F0126 


< 600 


0.90 


60 


45° 


<1.4 10^ 


<3.5 10^ 


2203 -F29 


< 1500 


0.67 


60 


45° 


<4.8 10^ 


<1.4 10^ 


2214 -M39 


800 


0.13 


30 


0° 


2.5 10^ 


2.6 10^ 


2235 -0301 


< 600 


1.55 


30 


0° 


<2.1 10^ 


<6.0 10^ 


2237 -0607 


< 400 


1.97 


30 


0° 


<1.5 10^ 


<7.4 10^ 


2331 -F0216 


< 1000 


0.87 


30 


45° 


<1.8 10^ 


<5.2 10^ 



TABLE 4 
Damped Lv-a systems 



system 


Zem 


'^abs 


[Fe/H] 


k/kc 


1020cm-2 


10-3 


Mq 


0000-26 


4.108 


3.3901 


> -2.44 


0.002 


25.7 


> 5 


> 2.4 10^ 


1202-0725 


4.69 


4.3829 


-2.23 


0.01 


4 


3.2 


0.48 10^ 


1331-M704 


2.08 


1.776 


-2.05 


0.015 


15.1 


22.7 


2.1 10^ 


2237-0607 


4.56 


4.08 


-2.18 


0.011 


3.3 


2.6 


0.44 10^ 



